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Abstract 

The scalar field degree of freedom in Einstein's plus Matter field equations is de- 
coupled for Bianchi type I and V general cosmological models. The source, apart 
from the minimally coupled scalar field with arbitrary potential V(3>), is provided 
by a perfect fluid obeying a general equation of state p = p{p). The resulting ODE 
is, by an appropriate choice of final time gauge affiliated to the scalar field, reduced 
to 1st order, and then the system is completely integrated for arbitrary choices of 
the potential and the equation of state. 
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1 Introduction 



In the past two decades there has been a growing interest in scalar field cosmological 
models primarily due to the prominent importance of scalar fields for inflationary sce- 
narios £Q. The importance of the coupling between a scalar field and the gravitational 
field has been further stressed by Madsen [2] , who has shown that it can have non-trivial 
consequences for the spontaneous breaking of gauge symmetries. A dynamical systems 
approach has been extensively used in the study of scalar field cosmologies and their 
asymptotic behavior See also [I] for a concise review. 

As far as exact solutions of scalar field cosmologies are concerned, Burd and Barrow 
[2] have studied homogeneous but anisotropic Bianchi models of types III and VI (as well 
as Kantowski-Sachs models) and have found exact solutions. Lidsey, and Aguirregabiria 
et al jBj have found exact solutions for Bianchi type I models. Feinstein and Ibanez [2j 
have found exact solutions for Bianchi models of type III and VI, Moss and Wright, 
Madsen, and Abreut, Crawford and Mimoso [E] have studied exact solutions in the 
setting of conformal scalar field cosmologies. Paul jU] has obtained exact solutions of 
a higher derivative theory in the presence of an interacting scalar field. The discovery 
of the BTZ black hole has motivated the study of analytic solutions in the context of 
scalar field cosmology in (2+1) dimensions [10] . while Russo ^1] has obtained the general 
solution for a scalar field cosmology in d dimensions with exponential potentials and a 
flat Robertson- Walker metric. An early work on exponential potentials has been done 
by Salopek and Bond (12] while recent treatments are given by Kehagias and Kofinas 
[T3] and Neupane [T3j . 

Cosmological models containing both a fluid and a scalar field have also been studied. 
Chimento and Jakubi [15J have given exact solutions of scalar field cosmologies with a 
perfect fluid and a viscous fluid respectively. Mendez ^H] has obtained an exact solution 
for the case of an imperfect fluid in a FRW spacetime. In the so called scaling scalar 
field cosmologies, the energy density due to the scalar field is proportional to the energy 
density of the perfect fluid. Thus a number of spatially flat, isotropic models in which 
the energy density of the scalar field is proportional to that of the perfect fluid have been 
investigated [Ej. Billyard, Coley and van den Hoogen [TH] have studied the stability of 
these scaling solutions within the class of spatially homogeneous cosmological models 
with a barotropic fluid matter content. Furthermore, they have studied the qualitative 
behavior of spatially homogeneous models with a barotropic fluid and a non-interacting 
scalar field with an exponential potential in the class of Bianchi type B models |19j . 
Saha [20] has obtained exact solutions for a Bianchi type I model with a perfect fluid 
and dark energy content, while Chimento and Cossarini have studied exact solutions in 
1+1 dimensions using an isotropic perfect fluid source [2*T] . 

In a significant paper, Hawkins and Lidsey [22], have shown that, for a flat FRW 
geometry, the dynamics of scalar field cosmologies with a perfect fluid matter content 
can be described by the non-linear, Ermakov-Pinney equation (which, in turn, leads to 
tantalizing analogies with the dynamics of Bose-Einstein condensates), while an early 
work of Barrow also deserves mention in this context. Exact solutions have been obtained 
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in this description [23]. Another kind of decoupling of the scalar field degree of freedom 
has been initiated in [21], where a Robertson- Walker background geometry minimally 
coupled to the scalar field has been investigated. By use of integrals of the motion and 
of the Klein-Gordon equation in the quadratic constraint equation, a single, higher-order 
non-linear differential equation for the scalar field was obtained. 

In the present work we generalize this decoupling for the case of general Bianchi type 
I and V geometries in the presence of a general perfect fluid source. The advantage is 
that the resulting ODE is fully integrated and this is achieved for arbitrary choices of the 
scalar field potential and the fluid equation of state. The paper is organized as follows. 
In section 2 the geometry and the matter description as well as the governing equations 
for the system are presented. In section 3 the decoupling of the scalar field degree of 
freedom is performed for Bianchi type I, the reduction of the resulting ODE is given and 
the system is completely integrated. The corresponding calculations for Bianchi type V 
are given in section 4. Finally, the conclusions and a discussion of the results obtained 
are presented in section 5. 

Throughout we use geometrized units, i.e. c = 8nG = 1, while has the signature 
(-,+,+,+). 



2 The governing equations and the matter content 

Our starting point is the line element for a spatially homogeneous geometry: 

ds 2 = [N a (t)N a (t) - N 2 (t)] dt 2 + 2N a (t)a°(x) dx l dt + 7a/3 (t)af (x)af (x) dx l dx j (2.1) 

where crf(x) are the basis one- forms satisfying ofj(x) — of^x) = 2Cg cr| \x)a] \x) , (with 
C^ 7 the structure constants of the Lie algebra of the corresponding three- 
dimensional isometry group of motion acting simply transitively on the spatial hyper- 
surfaces of simultaneity), N a (t) is the shift vector, N(t) is the lapse function, and 7 Q( g(t) 
is the scale factor matrix of the Bianchi type examined. 

As explained in [2H], there are special general coordinate transformations mixing 
space and time in the new space coordinates, whose effect on the line element (12.1)1 is 
described by 

%p(t) = K^u(t), 

N(t) = N(t) (2.2) 
N a (t) = Ai(t) [N p (t) + P0(t) lpP (t)} 

where A£(t), P p {t) satisfy 

K C M = k K c % ( 2 - 3 ) 

2P"C£,A£ = k% (2.4) 

the overdot denoting differentiation with respect to time. Due to (|2.3J) . A%(t) belongs 
to the automorphism group of the Bianchi type in question and transformations (|2.2j) 
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describe the gauge freedom of the emanating system of Einstein's field equations in 
the case of vacuum. So, the three arbitrary functions of time contained in A?(t) and 
P a (t) can be used to simplify the line element (j2.1J) and, thus, also the aforementioned 
equations. Automorphisms induced by general coordinate transformations have also 
been considered in j2E], while the rigid gauge symmetries have been analyzed in |27j . 
Time-dependent automorphisms, seen as tangent space transformations, have also been 
considered previously [?], (29, while the first ( known to us) reference of the relevance of 
Automorphisms to a systematic analysis of Bianchi Cosmologies goes back to 1962 [30J. 

The above result also holds true when the matter content is such that the linear 
constraints G® do not acquire an extra term Tf. For orthogonal perfect fluids and scalar 
fields depending only on time, Tf is indeed zero and therefore we can use the gauge 
freedom (j2.2H2.4|) to diagonalize 7 Q/ g(t) without loosing generality. Therefore, the scale 
factor matrix is taken to be 

la p{t) =diag(a 2 ,fo 2 ,c 2 ) (2.5) 

while for the shift vector we have N a (t) = 0. Our choice of time is specified by the 
gauge condition N = ^detj^, which is frequently used due to the simplification of the 
Einstein tensor and whose importance for the decoupling is essential. 

The matter content is a minimally coupled scalar field with an arbitrary potential 
V(§), thus having an energy-momentum tensor 

TjS = - i^O^^A + 2V(*)) (2.6) 

and a perfect fluid part 

Tjj$ = {p + p)u ll u v +pg va , (2.7) 

where m m is the unit four-velocity vector and a general equation of state p = p(p) is 
adopted. The governing Einstein's Field Equations are taken to be 

B!t - = Tj? (2.8) 

where = T$ + T$) , while the scalar field equation (Klein-Gordon) is 

1 



-g 



fy(V=0<r&$) - V($) = (2.9) 



the prime, from now onwards, denoting differentiation with respect to the argument. 

The "equation of motion" for the perfect fluid is the conservation of its energy- 
momentum tensor: 

= o (2.io) 

(The scalar field energy-momentum tensor is separately conserved by virtue of the Klein- 
Gordon equation). 
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3 Decoupling of the scalar degree and the solution 
space for Bianchi Type I 

The basis one-forms are crf{x) = 8f and, thus, with 1)2. 5 j) . zero shift and the chosen time 
gauge, the initial form of the metric is given by 
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(3.1) 
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Spatial homogeneity implies that $ = and p = p(t). The nonzero components of 
;he Eii 

a 2 b 2 c 2 ) are given by 



the Einstein tensor G^, and the energy- momentum tensors Tj; 1 , Tj 2 ^ (all multiplied by 
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Due to the equalities ()3.7|) and (|3.9j) . we can subtract the corresponding Einstein' 
equations = + Tj 2 ^, i.e. form the differences (all multiplied by a 2 b 2 c 2 



G\ 



G\ = T 2 (1)2 + 71 



,(2)2 



n 



(i)i 



T[ 2)1 and G\ 



T (l)3 + T 



(2)3 



and get the following equations involving only the scale factors 



b 2 a b 
b 2 a b 



T (l)l _ T (2) 



(3.10 
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^-£!-! + ! = (3.ii) 

ar c z a c 
These equations provide the following integrals of motion: 

b{t) = e xt a(t) (3.12) 



c (t) = e^a{t) (3.13) 

leaving us with just one undetermined scale factor. At this stage the Klein-Gordon 
equation (|2.9|) becomes 

~ V 'W 5^ = (3-14) 

Now, the conservation of the fluid energy-momentum tensor (continuity equation) gives 

3- + —-^ + A + ^ = (3.15) 

a p(p) + p 

It is straightforward to check that the time derivative of the quadratic constraint 
equation G® = Tq + Tq 2 ^ is identically satisfied by virtue of the remaining spatial 
equation G\ = T^ 1 + 1 , the Klein-Gordon equation and the continuity equation 
(solved for a(t), $(£), p(t)), as expected from the consistency between the aforementioned 
constraint with the spatial Einstein and the Matter equations. Therefore, the equations 
to be solved are the Klein-Gordon equation (|3.14|) . the continuity equation ()3.15|) and 
the constraint equation which, upon multiplication by — 2a 6 e 2( - x+ ^ t reads: 

6^-) +4(A + /i)- + 2A//-$ 2 -2e 2 ( A+ ^a 6 [^($) + p] =0 (3.16) 

In order to have a closed form for the integral of the continuity eq. ()3.15|) , it is convenient 
to use the parametrization 

pto = fg-p (»■"> 

(the prime denoting differentiation with respect to the argument) through the use of 
which one obtains the integral 

p = h[p a~ 3 e- ix+ ^} (3.18) 
with h being the inverse function to g, i.e. satisfying h(g(x)) = x. 

^From the Klein-Gordon equation ()3.14j) the scale factor and its logarithmic derivative 
is expressed in terms of the scalar field $(t) and its derivatives as well as V(<&): 

-2(A+M)t<bl 1/6 
7-t— (3-19) 
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a_l(S, (3.20) 



a 6 V$ V'{$) 

The promised decoupling of the scalar field dynamics from the geometry occurs upon 
inserting (j3.18H3.2T)|) into (jH.lfij) . The result is the following non-linear, third order ODE 
for $(£): 

2 



+12 vW) v ^ + 12 vW) h 



. > ~ 1/2 ' 



Po 



(3.21) 



At this stage, any solution to this equation determines, through ()3.18j) and ()3.19|) a cor- 
responding solution to the full Einstein plus Matter system whose entire space of solution 
is therefore attained from the solution space of ()3.21|) . Of course, the price paid, for the 
moment, is the non-linearity in the highest time derivative, which has also been raised 
to third order. Furthermore, normally one would expect that the arbitrary functions h, 
V need first to be specified before hoping to actually get a solution. Nevertheless, it is 
quite interesting that further reduction of the order of ()3.21|) . and subsequent complete 
integration of the whole system, is possible. To this end, we first observe that many 
terms in (|3.21j) . namely all the nontrivial except the 4th and 5th, are functions of the 
$ 

combination — — — -. We thus define: 

* - ~vW) (3 ' 22) 

and write ()3.21|) as 

2 



-4(A + ^) 2 + 12A/i-6/(x)-12 X %oX" 1/2 )=0 (3.23) 



where 



fix) = ^ 2 + 2xV($) (3.24) 

implicitly reflects the arbitrariness in choosing V A ($). Now ()3.23|) can be integrated and, 
by judicious choices for /, h, even give x(t) i n closed form. Suppose x(t) does indeed 
solve ()3.23|) . Then, multiplying ()3.22|) by 2$V( ( I ) ) and using ()3.24|) we get the first order 
linear differential equation for f(t): 

df(t) 



dt 



2 X (t)V(t) (3.25) 



and consequently $(t) is given by 

d$(t) 



(Jf ±Vf(x(t))-2 X (t)V(t) (3.26) 
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Finally, a change of time variable from t to x (defined in (|3.22jl ) permits the presentation 
of the entire space of solutions to the system under consideration in closed form. Indeed, 
considering V, f, $ as functions of X one gets from (j3.25H3.2()|) respectively 

T = 2V{x) (3 ' 27) 

^-X(t) = ±Vf(x)-2 X V( X ) (3.28) 

where the time derivative of x is given by the constraint ()3.23|) . The integration of these 
two equations is trivial, yielding /, $ as: 

/(*)=* + 2 [v(x)d X (3.29) 



* (X) = " * / \ V 4(A + - J£+ SuxHKXT*) dX (3M) 
The line element in the new time x reads: 

dx 2 



ds 1 



X [4(A + - 12A/x + 6/( X ) + l2 X h(p x- 1/2 )} 



(3.31) 



+ [e -2(A +M ) t (x) x ]l/3 dx 2 + [e (4A-2 M ) t (x) x ]l/3 dy 2 + [e (-2A+^)t(x) ^] 1/3 ^2 

with t(x) given by the integral form of ()3.23|) 

t( X )=± I 1 d X (3.32) 

while the density and pressure are given as: 

P(X) = h(poX~ l/2 ) P(X) = PoX- 1/2 h'(p oX - l/2 ) ~ h( PoX ~ l/2 ) (3.33) 

the prime denoting differentiation with respect to the argument. Quite independently 
of the way these solutions were obtained, one can straightforwardly check ( through,say, 
a symbolic computing facility) that they do satisfy all ten Einstein's Equations, the 
generalized K-G and the continuity equation. Furthermore, since no extra ansatz has 
been involved in the process of integration of the system (j2.8H2.10jl . equations (|3.29H3.33j) 
represent the full space of solutions to the Einstein plus Matter system considered. The 
functions V, h can be freely specified to obtain special case solutions. 
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4 Decoupling of the scalar degree and the solution 
space for the Bianchi Type V 

In this case the basis one-forms are a 1 = e~ x dy, a 2 = e~ x dz, a 3 = dx. The G\ = 
Einstein equation (due to spatial homogeneity $ = and p = p(t) and therefore 
there is no corresponding component of the matter tensor) implies c = ab. Thus the 
initial metric is taken as: 



9nv 
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The situation is similar to the Type I case, and thus forming the difference (multiplied 
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b 
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which can be integrated yielding 



b(t) = e xt a(t) 



(4.10) 
(4.11) 
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The Klein-Gordon equation (|2.9jl becomes 

e~ 3A '$(T) 

while the conservation of the fluid energy-momentum tensor (continuity equation) gives 

3 - + ^ + |A = (4.13) 

a p(p)+p 2 

Again, the only other equation to be solved is the constraint equation which reads 

6 ( -V + 6A- + A 2 - 6e 2A V - $ 2 - 2e 3Xt a 6 (V(<5>) + p) = (4.14) 
\a J a 

where the fourth term is the only non trivial difference from the corresponding equation 
()3.16|) and its presence is due to the non-vanishing curvature of the spatial slice. 

Integrating 1)4.13)1 (in the parameterization (J3.17)) ) and solving ()4.12j) for a(t), we 
obtain the following results for the matter density p(t), the scale factor a(t) and its 
logarithmic derivative 

p = h (p a~ 3 e-^ xt ) (4.15) 



-3At$ * ' " 



a 1 V"(<$>) 



(4.16) 

$ - 3A ) (4.17) 

3A 'a 6 ) 



a 6 V$ V'($) 

Use of these equations in the quadratic constraint equation (|4.14jl (multiplied by — 2e 
yields 

(I - m*V - 3A' - 36 (—L>r - 6i 2 + 12 * vm 



(4.18) 



+ 12 7777^ A> 



which, with the same definition of x, translates into 

3A 2 - 36 X 2/3 - 6/(x) - l2 X h(poX' 1/2 ) = (4.19) 



2 

X 



X 

By arguments completely analogous to the previous Type I case, the final form of the 
solution is, in this case: 

f(x) = a + 2 [ V( X )d X (4.20) 



®(x) 



K± 
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1 / f(x) - 2XV(X) 

X V 3A 2 + 36 X 2 / 3 + 6/( X ) + I2xh{p x- l/2 ) 



dx 



(4.21) 



dx 2 



X [3A 2 + 36 X 2 / 3 + 6/( X ) + 12 X ^(poX- 1/2 )] 



(4.22) 





Xv /3A 2 + 36x 2 / 3 + 6 f(x) + l2xh( PoX - 1/2 ) 



dx 



(4.23) 



while the density and pressure are given as: 



P(x) = HpoX 




P(x) = PoX 



/2 h'(poX- 1/2 ) ~ KpoX 



1/2 



(4.24) 



The remarks at the end of the previous section apply also here 

5 Discussion-Conclusions 

We have discussed the dynamics of a scalar field with an arbitrary potential, minimally 
coupled to a general (anisotropic) Bianchi Type I and V geometry, in the presence of a 
perfect fluid source obeying a general equation of state. In the case of vacuum, the rich 
structure of Outer Automorphisms for these two symmetry groups entails the existence of 
some integrals of motion: Indeed, consider the generators of the rigid symmetries ([2.3)1 . 
i.e. the vector fields in the space of dependant variables Xj = X p Ia j p p -^-^ where A^j 
satisfy \" p = X p I/3 + Aj C^ p . If one performs the appropriate change of variables 
la/3 z a/3 which brings one (or more, if possible) generator into each canonical form, 
say gj^-, then Einstein's field equations written in the new variables do not explicitly 
depend on z\\. Thus the system becomes of first order in the variable zu and therefore 
an integral of motion arises. 

The initial choice of the time gauge in which the lapse is equal to the determinant of 
the scale factor matrix has a twofold advantage: Firstly, it enables the corresponding 
integrals of motion (in the presence of the matter content chosen) to be revealed. Sec- 
ondly, it makes the Klein-Gordon equation purely algebraic in the scale factor variables. 
As a result, when the utilization of the integrals of motion has reduced the number 
of these variables to one, this equation gives this last scale factor as a function of the 
second derivative of the scalar field, the derivative of the potential with respect to the 
scalar field, and of time. Substitution of this form of the scale factor into the only equa- 
tion remaining to be solved, i.e. the quadratic constraint, results in a single ODE for 



10 



the scalar field (without any explicit time dependence!). The utilization of a final time 
gauge adapted to the scalar field enables the reduction of this equation to first order 
and subsequently leads to the complete integration of the entire system of Einstein plus 
Matter Field Equations. The description of the space of solutions contains, in an integral 
form, the arbitrary functions of the final time x-, V{x) and ^(poX -1 ^ 2 )- The presence 
of this twofold arbitrariness corresponds to the fact that we have not specified either 
the form of the potential (as a function of the scalar field) or the equation of state. It 
is evident that prescribing V(x) and h(poX 2 ) implicitly corresponds to a choice of 
potential form and equation of state. For example, the choice h = A(pQX~ 1 ^ 2 ) 1+ ' y gives 
through (|4.24jl p = jp, i.e. the barotropic equation of state. For the scalar field, if we 

take the particular case X — p, — a — /i(poX~^ 2 ) = (i- e -; flat Robertson- Walker with 

C 

no fluid) the choice V(x) = corresponds to the functional form V = ±Ce^( K ~®\ 

VX 

If, on the other hand, someone insists in prescribing V($>) and p(p) := p(h) then the 
situation must be dealt with in the following manner: As far as the density is concerned, 
the equation giving the pressure becomes the holonomic, first order differential equation 
p(h) = wh'(w) — h(w) with w standing for poX™ 1 ^ 2 which can be straightforwardly inte- 
grated. For example, a barotropic equation p := 7p gives p = A(pqx~ 1 ^ 2 ) 1+ " 1 ■ As far as 
the matter field is concerned, the situation is somewhat more complicated as the results 
of choosing a particular form V($) are influenced by the choice of the density h. As an 
example consider again the case A = p = a = ^(poX™ 1 ^ 2 ) = and an arbitrary potential 
form V(<&). Relaxing for a moment ([3.30)1 . we get the following Klein-Gordon equation 



V($) 6 X 2 $ /2 (x) 

which may be difficult to solve depending on the choice of V(<&). This is the price paid 
for insisting on prescribing V(&) and not V(x) i n which case the solution would be 
given by ()3 .29113. 33*j) or (|4. 2 0114.24^) correspondingly. The particular choice V = Vo e_A * 
corresponds to the case considered in [TTj and the equation above can be dealt with by 
choosing a new time x = e 3r/ ' 2 , in which the equation becomes the following first order 

m the derivative u = — — : 

dr 

48^ - 96cj 3 + 24Aw 2 - 9A = 
dr 

Finally, there are two cases which, at first sight, need to be separately examined. The 
first concerns the case of a constant potential V($) = Vq, i.e. a cosmological constant 
term. Then, V'($) = and the definition of the final time x m equation ()3.22|) seems 
to be precarious. However, the Klein-Gordon eq. 1)3.14)1 then implies that also $ = 
and, surprisingly enough, x does exist. Indeed, substituting l / ($) = Vq in the solutions 
or (I4.20H4.241) one can see that the solution is still valid. This holds true even 



$ 1 

for Vo = 0. The second case arises when the ratio - = — -, for, then, the change 

1/ '(<£>) & 
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of time from t to x is not valid. In this case, the remaining scale factor is given by 
a(t) = ce ^~ in the Type I case or a(t) = ce~% in the Type V case. The continuity 
equation ((jHHSJ) or (|4.12jl respectively) implies that the pressure and density are con- 
stants, say, p and p . The remaining equations dictate $(t) = At + B and V(<&) = V 
and there exists a relation between c,po,Po,A and B due to the quadratic constraint. 

The general properties of the space of solutions found for both type I and type V cases, 
e.g. isotropization, attractors and self-similarity, and quintessence, will be examined in 
a forthcoming paper, in which also particular cases of high interest are to be explicitly 
elaborated. Possible applications of the method exhibited here can be the cases of D + 1 
spatially homogeneous spacetimes, in which the Outer Automorphisms are rich enough 
to provide sufficient integrals of motion. Then the method will be applicable if the matter 
content is such that the integrals of motion persist. 
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